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‘ballistic  pressure  - "single  rc  md  parametric  resonance". 

(11)  the  periodic  applications  of  ballistic  pressure  such  as 
encountered  In  an  automatic  weapon  - "multiple  round 
parametric  resonance". 

Results  show  that  ballistic  cycles  currently  employed  In  the  60mm 
MCAAAC  semi  automatic  cannon  are  not  likely  to  excite  single  round 
resonance.  Unusually  brief  cycles,  however,  are  shown  to  be  capable 
of  producing  resonance  amplifications  of  three  orders  of  magnitude  In 
less  than  twenty  cycles  of  axial  vibration.  By  proper  design  of  the 
pressure  cycle  and/or  the  fundamental  axial  frequency  of  the  tube,  this 
type  of  resonance  is  rather  easily  avoided.  ^ 

Further  results  show  that  for  the  20mm  Ml 39  machine  gun,  ampli- 


fications in  excess  of  fifty  can  be  reached  in  under  five  seconds  of 
continuous  firing.  A special  application  of  the  work  of  Krajcinovic 
and  Herrmann  leads  to  a set  of  instability  contours  from  which  the 
growth  (characteristic)  exponent  can  be  determined  as  a function  of  the 
ratio  of  natural  and  excitation  frequencies  and  the  product  of  the 
ballistic  impulse  and  the  tube  slenderness  ratio.  Control  or 
elimination  of  multi-round  resonance  can  be  maintained  either  through 
control  of  the  initial  conditions  or  by  designing  for  mismatch 
between  the  transverse  frequencies  and  integral  mulitples  of  one-half 
the  excitation  frequency,  i.e.,  the  firing  rate. 
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INTRODUCTION  AND  BACKGROUND 


Forced  vibrations  of  undamped  linear  systems  are  characterized 
by  the  differential  equation: 

x + w?x  = f ; tuJ  = k/m  (1 ) 

where  m and  ,k  are  the  inertl;  and  stiffness  parameters  of  the  system 
and  x represents  the  system  displacement.  Conventionally  m and  k are 
constants  and  f Is  a time-variant  force  which  causes  resonance  if  it 
contains  a component  having  the  system  period  2ti  / /k/m  . The  resonant 
term  is  linear  in  the  time  variable  t and  is  a particular  solution 
of  (1). 

If  ai  Is  time  dependent  the  solution  of  equation  (1)  is  much  more 
complicated  and  in  most  cases  has  only  been  achieved  through  approximate 
methods  or  numerical  quadrature.  An  important  subclass  of  problems 
exists,  however,  for  which  a good  deal  of  theoretical  progress  has  been 
made.  These  are  problems  In  which  the  variation  of  w is  periodic  and 
f is  identically  zero.  Such  cases  are  representpd  by  the  homogeneous 
linear  differential  equation: 

x + v2(t)x  = 0 (2) 

where  w2  = i'2(l  - e<J>(t)) 
o 

and  $ is  periodic  in  time.  Since  this  equation  is  homogeneous  it  admits 
a general  solution  of  the  form: 

x * Ax-j  (t)  + Bx2( t ) (2a) 

where  A and  B depend  only  on  the  initial  conditions  of  the  problem. 
Floquet's  theorem^  allows  for  two  solutions  of  the  form: 

^Morse,  P.  M.  arid  Feshbach,  H.,  Methods  of  Theoretical  Physics,  McGraw 
Hill,  1953,  557.  - - 
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(3) 


x ^ ( t ) = <j»1eY  t 
x2(t)  = 4>2e'Yt 

Historically,  equation  (?)  is  known  as  Hill's  Equation  and  (3)  are 
its  Floquet  solutions.  (See  Appendix  A.)  The  4>^(t)  have  the  same 
periodicity  as  the  'excitation'  function  $ (t).  Thus  if  y has  a non- 
zero real  part,  one  of  these  solutions  is  unstable  and  the  general 
solution  exhibits  exponential  growth  provided  that  the  Initial 

conditions  are  not  those  which  would  cause  the  corresponding  coefficient 

2 

of  the  growth  term  to  vanish.  Theoretically  it  has  been  shown  that 
unstable  solutions  car  result  whenever  the  ratio  of  a system  natural 
frequency  to  the  frequency  of  excitation  takes  on  values  in  the 
neighborhoods  of  Integral  multiples  of  one-half  (cf  rig.  1).  Thus  the 
primary  Instability,  for  example,  will  be  encountered  when  the 
excitation  frequency  approaches  twice  a natural  frequency  of  the 
system.  We  therefore  have  three  fundamental  differences  between 
conventional  forced-resonance  and  that  induced  parametrically: 

(a)  Forced  resonance  is  Independent  of  the  initial  conditions 
whereas  parametric  resonance  is  not.  Given  a force  component 
operating  at  a natural  frequency  of  a system,  resonance  must  occur 
whereas  a parameter  'stiffness,  mass)  varying  periodically  at  2/n 
times  the  system  frequency  (n  an  integer)  need  not  produce  resonance 

2Bolotin,  V.  V.,  The  Dynamic  Stability  of  Elastic  Systems,  Holden-Day, 
1%4,  22-23. 


If  the  Initial  conditions  can  be  controlled.  This  Is  especially  signif- 
icant when  the  short  term  response  of  a system  is  of  Interest  for  there 
is  a wide  choice  of  initial  conditions  for  which  the  early  response  will 
have  decreasing  amplitude. 

(b)  Forced  resonance  consists  of  oscillations  whose  amplitude 
increases  linearly  with  time  whereas  parametric  resonance  produces 
exponential  growth. 

(c)  Forced  resonance  occurs  if  and  only  if  the  forcing  frequency 
exactly  equals  a natural  frequency  of  the  system.  In  contrast  parametric 
resonance  can  occur  whenever  an  Integral  multiple  of  the  excitation 
frequency  approaches  twice  the  value  of  a natural  frequency.  That  Is, 
parametric  resonance  - unlike  forced  resonance  - Is  not  a singular 
phenomenon  but  occurs  throughout  the  neighborhood  regions  of  a countable 
infinity  of  critical  frequency  ratios.  It  is  therefore  a regional 
phenomenon.  An  Infinity  of  unstable  regions  exist,  the  most  important 

of  which  is  the  primary  region  of  instability. 


Effect  of  Linear  Dampinc 


The  addition  of  the  linear  damping  term,  2cx,  Into  equation  (2) 


creates  no  complication  since  a transformation  x * vw  can  always  be 

- /cdt 

found  (even  when  c is  time  dependent)  - such  that  w = e and  v solves 


A 


the  differential  equation: 

v - [c2  + c - u)2(l  - e$(t))]v  a 0 

• • 

which,  if  c = 0 , has  the  form  of  equation  (2)  and  hence  has  solutions 


Thus  e^ft-Jcdt]  [yt+/cdt] 

x = vw  = A<^  (t)  + B<J>2(t) 


■i  ■ 


i.e.,  the  inclusion  of  a linear  damping  term  results  in  a simple 


Two  Examples  of  HI  1 1 ' s Fquation 
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several  examples  of  hITT's  equation  art-  given  in  the  literature  ’ 
The  one  most  often  cited  pertains  to  a beam  column  subjected  to  a 
periodically  varying  axial  compressive  load  P(t)  such  os  depicted  in 
Figure  a.  below: 


Figure  a. 


ilassical  Beam  Problem  Governed  by  Hill's  Equation 


The  govern.nq  differential  equation  from  Fuler-Bernoul li  beam 
theory  is: 


El 


iX 


Pit)  ^ « 0 

ix*  Tt* 


(4) 


As  Fraicinovic  and  Hermann3  have  oointed  out,  an  attempt  to 
separate  variables  through  the  substitution  u(x,t)  = X(t)f(t)  will 
result  in  the  ordinary  differentia'  equation: 


-Bolotin,  V.  V.,  The  Dynamic  Stability  of  Elastic  Systems,  Holden-Day, 
1964,  22-23.  ' ~~ 

3Den  Hartog,  J.  P.,  Mechanical  Vibrations,  McGraw  Hill,  1940,  378. 

5«ra  jc  movie  , r.  ►’ . and  Hermann,  G.  , S tab i 1 i ty  of  Straight  Bars 
5 j b je cted  t o_ Re oea ted  J mpu Is i ve  Comp res s Ton , AIAA  Journo T OcT  68 , 

2u2S-2C‘7.'' 


Since  the  right  hand  side  of  this  equation  is  independent  of  x: 


HA  = const  = a 
P X 


and 

X " 

— = const  = aB 

PX 

Thus 

f + cx[l  - BP(t)]f  = 0 

(5a) 

and 

ElXiV  - apX  = 0 

(5b) 

and 

EIX1V  + 1 X"  . 0 

(5c) 

While  eauation  (5a)  is  the  desired  Hill's  equation  for  the  system, 
the  separation  of  variables  approach  is  only  valid  when  X(x)  also 
satisfies  equations  (5b)  and  (5c);  i.e.,  the  modes  of  free  vibration 
which  solve  (5b),  must  be  identical  to  the  buckling  modes  which  solve 
(5c).  The  case  depicted  in  figure  a.  - a hinged-hinged  support 
system  - does  in  fact  satisfy  both  of  these  conditions.  In  most  cases, 
however,  the  boundary  conditions  lead  to  modes  which  do  not  satisfy 
both  (5b)  and  (5c).  In  such  cases  an  approximate  Hill's  equation  can 
be  obtained  through  a variacional  procedure,  such  as  that  due  to 
Galerkin.  In  either  case,  therefore,  the  problem  Is  reduced  to  the 
analysis  of  equation  (5a)  where  a and  B derive  from  an  analysis  which 
is  either  exact  or  approximate. 

Another  example  - one  which  more  directly  leads  to  Hill's 
equation  - considers  a system  with  time-variant  inertia,  such  as  a 
child  pumping  a swing  (figure  b).  Essentially  a concentrated  mass  is 
raised  and  lowered  periodically  along  a relatively  massless  rod 


i 


f 


(or  chain,  etc.)  which  pivots  at  0.  The  rotational  Inertia  of  the  j 

I 

pendulum  thus  varies  periodically  In  time.  The  example  is  one  where 
the  path  taken  by  the  mass  through  the  gravitational  field  results  in 
a net  amount  of  work  per  cycle  of  the  motion. 

The  equation  of  motion: 


mg 


-S,(t)mg  sine  = ^ (mJt2(t)§) 


. , dit  de  , d*6 

' 2m  eft  dt  dt2  (6) 


Approximating:  sin  e 

(U  de 
dt  dt 


e. 


+ 2^  — +q6  = 0 
dt2  *+ 


ri cure  h.  - Swing  Problem  Leading  to 
Hill ' s Equation 


Defining  Jt  = lQ+h(t),  t - wt,  and  w2  = g/Jt0  leads  to  the  nondlmens^onal 
equivalent  of  equation  (6): 

0 + Mil)  0 + lhe  + e = O ; (')  = d/dx 

Iq  Iq 

If  h(x)  = ZQe  cos  2x  and  only  first  order  terms  in  e are  retained: 

0 - 4e  sin  2x0  + (1  - e cos  2x)e  = 0 (7) 
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I \J.  j 


Through  the  afore-mentioned  transformation  6 * vw,  the  equation  for 
v and  w are: 

v + (1  + 3e  cos  2t)v  = 0 (Hill's  Equation) 


and 


w . e-e  cos  2x 


A grjat  amount  of  consideration  has  been  given6  to  cases  where  the 
'characteristic  exponent'  y,  is  purely  imaginary  and  the  periodic 
excitation  <i>  (t)  is  sinusoidal.  Equation  (2)  then  reads: 

x + (a  - 2q  cos  2t)x  = 0 (8) 

This  equation  - a special  case  of  Hill's  equation  - is  called 
Mathieu's  Equation  (canonical  form).  As  with  any  Hill's  equation, 
Mathieu's  equation  yields  periodic  solutions  (called  Mathieu  functions) 
corresponding  to  purely  imaginary,  rational  values  of  the  characteristic 
exponent  y.  With  a view  toward  special  armament  applications,  however, 
this  report  will  deal  only  with  the  unstable  solutions  of  Hill's 
equation,  i.e.,  those  cases  in  which  y is  real. 

PARAMETRIC  EXCITATION  - ARMAMENT 

There  are  at  least  two  possible  sources  of  parametric  excitation 
in  gun  tubes  - that  is,  two  ways  in  which  periodic  coefficients  can  be 
introduced  into  the  beam  equations  of  motion.  The  most  obvious  can  be 
called  'multiple  round  excitation'  and  derives  from  the  periodicity 
present  in  automatic  weapons  in  which  several  time-variant  forces 
operate  at  the  firing  rate  of  the  weapon.  A reasonably  comprehensive 

6McLachlan,  N.  W. , Theory  and  Application  of  Mathieu  Functions, 

Oxford  Clarendon  Press,  1947. 
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differential  equation  including  these  forces  was  derived  in  a previous 
report^.  Figure  2 shows  a cantilevered  beam  model  of  a gun  tube 
acted  upon  by  several  curvature- induced  loads  (constant  projectile 
velocity  is  assumed  for  simplicity).  In  general  one  observes 
several  time  dependent  coefficients  multiplying  the  various  displace- 
ment derivatives.  In  automatic  weapons,  these  coefficients  are 
reproduced  periodically  according  to  the  firing  rate  and  will  appear 
in  the  Hill's  equation  obtained  upon  integration  of  the  space  variable. 
In  this  report  the  effects  of  only  one  such  term  will  be  investigated  - 
namely,  that  corresponding  to  the  periodic  ballistic  pressure  applied 
axially  at  the  breech. 

A second  and  less  obvious  cause  of  parametric  excitation  derives 
from  the  coupling  between  axial  and  transverse  tube  vibratio.is.  The 
simplest  equation  incorporating  the  necessary  nonlinear  coupling  terms 

O 

was  derived  by  Mclvor  and  Bernard0  In  1973.  Essentially  the  idea  Is 
that  a single  impulsively  applied  load  will  set  a column  ringing  with 
free  axial  vibrations.  Nonlinear  terms  - oscillating  at  the  frequency 
of  these  vibrations,  couple  with  the  transverse  displacement  variables 
chrough  the  stiffness  coefficients.  We  can  call  this  'single  round 
excitation'.  Thus  kinetic  energy  from  the  axial  vibrations  can  feed 
transverse  modes  and  lead  to  parametric  resonance.  The  governing 

7Simkins,  T.,  Pflegl,  G. , Scanlon,  R. , Dynamic  Response  of  the  Ml  1 3 
Gun  Tube  to  Travelling  Ballistic  Pressure  and  Data  Smoothing  as  Applied 
to  XM150  Acceleration  Data,  WVT-TR-75015^ 

8McIvor,  J.  K. , and  Bernard,  0.  E.,  The  Dynamic  Response  of  Columns 
Under  Short  Duration  Axial  Loads,  Trans  ASKiE,  September  1973,  688. 
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differential  equations  were  solved  by  the  authors  through  a Galerkln 
procedure  for  the  special  case  of  a simply  supported  beam  subjected  to 
an  axial  end  load  of  short  duration.  [It  should  be  noted,  however, 
that  there  Is  no  guarantee  that  the  variational  quantity  employed  will 
indeed  admit  an  extremum  when  the  associated  differential  operator  is 
nonlinear.]  Since  a good  deal  of  energy  is  apt  to  be  cops;. ,.ed  in 
rigid  body  recoil  in  armament  applications,  fixed  supports  are  to  be 
avoided.  Consequently  a tube  (beam)  cantilevered  from  end  supports 
which  allow  axial  movement  (Figure  5)  was  chosen  as  the  subject  of 
analysis  for  this  report.  (Relative  motion  of  the  supoort  is  ignored.) 
Evidence  of  Parametric  Resonance  in  Gun  Tubes 

In  order  to  minimize  shot  dispersion  in  automatic  weapons  the 

Q 

current  design  handbook  dealing  with  gun  tube  design  advises  that 
the  ratio  of  the  fundamental  transverse  frequency  of  the  tube  to  the 
firing  rate  be  kept  greater  than  3.5.  The  basis  for  this  value  is  a 
plot  of  shot  dispersion  vs.  frequency  ratio  R^  appearing  in  the  handbook 
and  reproduced  as  figure  3 of  this  report.  Referring  to  this  figure 
three  very  prominent  maxima  are  observed  at  successive  integral  values 
of  Rf  = 1,2  and  3.  The  reference  cited  in  connection  with  this  plot 
is  a 1955  report  by  Wente,  Shoenberger  and  Quinn  o^  Purdue  University^. 
Their  results,  shown  in  figute  4,  are  in  marked  contrast  to  those  of 
figure  3,  however.  Absent  is  the  maximum  at  Rf  = 3.0  shown  in  figure  3. 

^AMCP  706-252,  Engineering  Design  Handbook,  Gun  Series,  Gun  Tubes, 
February  1964. 

^Wente,  B.  E.,  Schoenberger,  R.  1.,  and  Quinn,  B.  E.,  An  Investigation 
of  the  Effort  of  the  National  Frequency  of  Vibration  of  the  Barrel 
Upon  the  Dispersion  of  an  Automatic  Weapon,  Purdue  U.,  1955,  AD64132. 
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It  Is  also  noted  that  figure  4 contains  no  information  below  the  value 
Rf  = 1.0.  Thus  the  only  features  common  to  both  figures  are  apparently 
the  maxima  at  - 1.0  and  = 2.0.  Accepting  these  maxima  as  the 
only  credible  information  to  be  gleaned  from  the  reference  publications 
one  searches  for  an  explanation  as  to  their  cause.  While  the  maximum 
at  Rf  = 1.0  may  be  attributed  either  to  parametric  or  to  ordinary 
(forced)  resonance,  that  at  = 2.0  cannot  be  due  to  ordinary 
resonance  and  may  be  evidence  of  parametric  resonance  - which,  as 
previously  discussed,  can  be  expected  to  occur  near  nominal  values  of 

Rf  - 1/2,  1,  3/2, n/2, Though  parametric  resonance  should  also 

produce  a dispersion  maximum  at  R^  = 1.5  in  Wente's  plot.  It  may  be 
that  it  has  been  missed  due  to  the  paucity  or  data  points. 

Equations  of  Motion 

The  model  chosen  to  represent  armament  applications  is  shown  in 
figure  5.  The  equations  of  motion  which  include  coupling  between 

O 

transverse  and  axial  displacements  are  those  of  Mclvor  and  Bernard  . 
Eliminating  their  dissipation  parameter  for  simplicity,  these  are: 


U-[u'  + 1 v'2]'  =0 
2 

(9a) 

v + azvi v - (u V ) ' = 0 

(9b) 

Boundary  Conditions 
P(t)/EA  = ea(0,t) 


0 = ea(l  ,t) 


L3s  2 Us>  Jo,T 
+ 1 (IVv  2-, 

Us  2 Us;  J1,T 


a 

b 


8McIvor,  J.  K. , and  Bernard,  J.  E.,  The  Dynamic  Response  of  Columns 
Under  Short  Duration  Axial  Loads,  Trans  ASME,  September  1 $73,  688. 
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0 ■ v(0,t)  C 
0 = v" (1  ,t)  d 
0 = v'(C>x)  e 
n = v"  ' (0,t)  f 


where:  u - n/L  = dimensionless  axial  displacement 

v - £/L  = dimensionless  transverse  displacement 
s = x/L  = dimensionless  coordinate 

a2  * I/AL2  = square  of  the  reciprocal  of  the  slenderness  ratio 
t * at/L  = dimensionless  time  variable 
a = (EA/p)1^2  = extensional  wave  speed 
E » Young's  Modulus  of  Elasticity 
A » Beam  Cross  sectional  area 
p = mass  per  unit  length 
I = area  moment  of  inertia 

P(t)  * end  loading,  a ballistic  pressure  function  of 
duration  tq 

e = axial  strain  Including  lowest  ordered  nonlinearity 
d 

Multiplying  (9a)  and  (9b)  by  <5u  and  6\>  respectively  and  Integrating 
over  the  length  of  the  beam: 

/' (U6u  - [u’  + i v'2]'5u)ds  = 0 (10a) 

0 

/ (V6v  - [ u ‘ v 1 3 1 5 v + a2v^v6v)ds  = 0 (10b) 

0 
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Using  boundary  conditions  a and  b together  with  an  integration 
by  parts,  (10a)  becomes: 

/'{ USu  + [u*  + 7 v,2]6u'}ds  + 6u(0,t)  = 0 (11a) 

0 1 tH 

Similarly,  boundary  conditions  b through  f applied  to  (10b)  yield: 

/ (V6v  + u'v'6v'  + a2v"6v")ds  + \ v'  (1  ,t)3i$v(1  ,t)  * 0 (lib) 

0 2 

The  last  term  in  (lib)  Is  an  order  higher  than  those  retained  and 
is  therefore  ignored.  Except  for  the  boundary  terms,  equations  (11) 
are  identical  to  those  obtained  by  Mclvor  and  Bernard  for  the  case  of 
a simply  supported  beam  under  end  loading.  Our  boundary  conditions, 
however,  dictate  a cnoice  of  completely  different  approximating 
functions  in  the  Galerkln  approach  to  solution.  Since  both  ends  of 
the  beam  in  our  problem  are  free  to  move  axially,  we  must  choose 
functions  which  do  not  constrain  the  function  u at  the  end  points  (0,1), 
i.e.,  there  are  no  geometric  constraints  such  as  are  present  v>  the 
problem  solved  by  Mclvor  and  Bernard.  A set  of  functions  which  appear 
to  satisfy  these  requirements: 

u = ^(t)  cos  jtts  ; j = 0,1,2,...  (12) 

j J 

For  the  transverse  motion,  the  eigenfunctions  of  a cantilevered  beam 
satisfy  the  conditions  c thru  f.  Hence: 

v = ITm(T)Wl71(s)  ; m=  1,2,3,...  (13) 
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♦ •.  • • • 


where 


W (s)  = coshes  - cosB^s  - ~ h&m  — cos^.  (slnh^s  - slnt^s) 
m ^ ™ slnhf^  + sln^  ™ ^ 

and  the  3^  have  values  1.875,  4.694,  7.855,  etc.;  further  values  can 

be  found  In  any  standard  vibrations  textbook  (cf  ref  11). 

Substitution  of  (12)  and  (13)  Into  equations  pi)  and  making  use 

of  the  orthogonality  of  the  trial  functions  and  the  Independence  of  the 

variational  quantities  6qj  and  6Tm  leads  to  the  following  sets  of 

nonll nearly  coupled,  ordinary  differential  equations: 


• • 
o 

II 

■n 

q0  * -P(t)/EA 

(14) 

j = 1,2,... 

Qj  ♦ (d*)2qj  + l AjmnTmTn  = -2P(t)/EA 
m,n 

(15a) 

m = 1,2,... 

Tm  * *VTm  Ajmn<jTn  = 0 

O.n  ^ 

(15b) 

where  the  A 

/ v'26u'ds 
0 


jmn  are  defined  ^rom  the  Integration: 


S T A-  T T 6q 
L jmn  m n M 

j.m.n 


j 


Table  I gives  the  values  of  these  coupling  coefficients  through 
A666*  Ration  (14)  Is  decoupled  from  the  others  and  being  repre- 
sentative of  rigid  body  motion  Is  of  no  further  interest.  Equations 

(15)  can  be  solved  numerically  by  any  of  several  numerical  Integration 
1 2 

programs  once  a has  been  specified  along  with  the  nondimensional 
load  function  P(t)/EA. 


^Nowackl,  W. , Dynamics  of  Elastic  Systems,  Chapman  and  Hall,  1963,  122. 

^Ralston  and  Wilt,  Mathematical  Methods  for  Digital  Computations, 

Wiley,  1960,  95-10$: 
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RESONANCE 


During  parametric  resonance  certain  transverse  coordinates  Tm(x) 
become  much  larger  than  the  rest.  Thus  the  quadratic  coupling  terms  can 
be  ignored  except  those  expected  to  resonate.  In  single  round  resonance 
it  is  expected  that  the  response  to  P(t)  will  mainly  occur  in  the  q^t) 
variable  (i.e.,  the  fundamental  axial  model  and  therefore  only  one 
such  term  need  be  considered.  If  m = M represents  the  particular 
transverse  mode  such  that  affy2  = tt/2,  primary  parametric  resonance  will 
result.  That  is,  the  natural  frequency  afy|2  Is  half  the  frequency  of 
the  exciting  variable  q^  = Asimrx,  t _>  xQ,  where  xQ  is  the  duration 
of  the  load  pulse  P(t). 

for  the  study  of  single  round  parametric  resonance  therefore, 
equations  (15)  reduce  to: 

+ Tr2q1  * -2P(t)/EA  = -2P*(t)  ; t i xQ  (16a) 

!h  * <“!V  * ■ 0 (16b) 

In  multiple  round  resonance,  transient  axial  vibrations  are 
ignored  and  the  qj(x)  are  assumed  to  follow  periodic  applications  of 
the  load  function  quasi-statically.  In  this  case,  if  the  firing  rate 
is  twice  the  toth  transverse  frequency,  parametric  resonance  will 
result.  Since  aV]_  of  the  q.(x)  are  periodic  according  to  the  firing 

J 

rate,  it  is  not  acceptable  to  retain  only  the  term  in  q^  as  in  the  case 
of  single  round  resonance.  However,  the  amplitudes  of  these  quasi - 
static  responses  attenuate  as  1 / j ^ (neglecting  the  quadratic  term  In 
15a)  and  only  a few  need  be  retained  for  accuracy.  In  place  of 
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equations  (15)  we  therefore  have: 


Qj  ♦ - -2P(t)/EA  r -2P*(t) 

(17a) 

rM  + (aV  + ?AJMM<,J(t))TM  * 0 

(17b) 

where  agM2  * (weapon  firing  rate)/2 

2tt 

The  periodic  character  of  the  q.  In  either  (16b)  or  (17b) 

J 

qualifies  these  as  Hill's  equations.  As  earlier  Indicated,  the 
general  solution  to  Hills'  equation  can  be  written  as: 

Tm(t)  = a<|>1(t)e>T  + b$2(T)e_YT  (18) 

where  the  are  periodic  functions  having  the  period  of  qj. 

It  should  be  mentioned  In  passing  that  when  y Is  real,  a plot  of 
the  solution  (18)  will  always  show  oscillations  near  the  natural 
frequency  aBM2  even  though  the  ^ have  the  same  period  as  q,  i.e., 
approximately  half  the  natural  period  . That  Is,  the  periodicity 
is  not  representative  of  the  oscillatory  appearance  of  the  response. 

It  can  be  shown  (see  Appendix  A)  that  once  a solution  for  Hill's 
equation  is  obtained  over  one  period  of  the  excitation,  it  can  be 
extended  analytically  for  all  time  by  means  of  Floquet's  theorem. 
Further,  in  the  case  of  Mathieu's  equation,  analytical  methods  have 
been  developed  leading  to  the  direct  determination  of  y,  the  character- 
istic exponent.  A detailed  series  of  curves  for  this  purpose  were 

■ o 

developed  by  S.  J.  Zaroodny10  in  1955. 

^3Zaroodny,  S.  J.,  An  Elementary  Review  of  the  Mathieu-HHI  Equation  of 
a Real  Variable  Based  on  Numerical  Solutions,  Ballistic  Research 
Laboratory  Memo.  Report  878,  Aberdeen  Proving  Ground,  MD,  1955. 
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EXAMPLE  I - THE  MCAAAC  GUN  TUBE 


It  is  of  Interest  to  assess  the  likelihood  of  encountering  single 
round  parametric  resonance  in  the  particularly  long  and  slender  tube 
planned  for  the  Medium  Caliber  Anti-Armor  Automatic  Cannon*  currently 
in  the  design  stage.  To  this  end  equations  ( l 6 ) will  be  employed  after 
first  establishing  the  magnitude  of  P*(t)  and  the  mode  M in  which 
resonance  might  be  expected. 

Equation  (16b)  always  possess  the  (trivial)  solution  * 0. 

This  solution  only  applies  when  the  'initial'  conditions  of  displacement 
and  velocity  following  the  application  of  the  load  P(t)  are  Identically 
zero  in  which  case  the  response  will  be  null  no  matter  how  intense 
the  'excitation'  Qj(t).  This  is  of  little  concern  in  armament 
applications  since  a good  deal  of  transverse  motion  is  certain  to  be 
excited  by  the  firing  of  a **ound.  For  example  the  recoil  of  a slightly 
curved  gun  tube  or  the  motion  of  the  projectile  therein  will  always 
excite  some  non-zero  'initial'  motion.  The  axial  vibration  q-j(T)  In 
response  to  the  ballistic  pressure  pulse  will  generally  result  in 
ampl ifica.ion  of  .hese  initial  motions  by  its  appearance  in  equation  (16) 
if  the  parameter  . happens  to  be  'tuned'  tor  parametric  resonance. 

According  to  results  from  the  latest  finite  element  (NASTRAN)  model 
of  the  MCAAAC  tube  (see  Tabic-  II),  the  fundamental  axial  frequency  is 
very  close  tc  being  twice  the  frequency  of  the  fifth  transverse  mode  - 

* The  MCAAAC  concept  plans  for  a two  or  three  round 'burst'  and  is  there- 
fire  not  an  automatic  weapon  in  the  same  sense  as  a conventional  machine 
fun.  Thus  multi -round  resonance  Is  thought  not  to  apply  and  the 
investigation  is  confined  to  that  of  resonance  which  might  be  induced 
from  the  firing  of  a single  round. 
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making  primary  parametric  resonance  of  this  mode  a possibility  If 
enough  axial  excitation  Is  produced  from  the  application  of  ballistic 
pressure  to  the  breech.  The  parameters  for  a uniform  beam  model  of 
this  tube  were  chosen,  therefore,  so  that  the  fundamental  axial 
frequency  of  580  hz  corresponds  to  ir  In  equation  (16a)  and  exactly 
half  this  value  corresponds  to  01852  of  equation  (16b).  A summary  of  the 
pertinent  parameters  Implied  by  these  assumptions  appears  In  Table  III. 
The  load  function  Is  approximated  by  a haverslne  shape.  Thus  (16a) 
reads : 


,» 

+ 7r2q.,  - -PJ(1  - cos  ~)  ; t < tq 

The  response  following  the  termination  of  ballistic  pressure  Is 
sinusoidal  with  amplitude: 


q-|  (max)  = 


2p*  /Sln7rTo/2v 

# _ j T/>  . . • 


Assuming  negll gable  response  from  the  other  axial  modes,  equation 
(16)  becomes  (for  suitable  choice  of  'initial'  time  zero): 

Tg  + a2854(l  + e cos  TrOTg  = 0 

where  e = A^q^  (max)/a26g4 


_2 

Using  values  from  Tables  I and  III,  e is  evaluated  at  2.07  x 10  . 
The  solution  for  T5  is  given  by  expression  (18).  In  view  of  the  small- 
ness of  e and  the  precise  state  of  tuning  assumed,  a very  good 
approximation  for  y can  be  obtained  by  the  method  of  strained  para- 
meters1*. The  result  Is: 

Y = tte/8 


14Nayfeh,  A,  H.,  Perturbation  Methods,  John  Wiley,  1973,  63. 
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In  real  time  the  magnitude  of  this  exponent  (',ee  Table  III)  becomes: 

VS.  x 1160  = 9.44 
8 

In  general  the  coefficient  of  the  growth  term  in  (18)  depends  on 
both  the  initial  displacement  and  velocity  of  the  Mth  transverse  mode 
at  'time  zero'.  It  is  this  quantity  which  is  amplified,  the  remaining 
term  of  (18)  becoming  less  important  as  time  moves  on.  As  a specific 
example  it  is  notable  that  for  certain  initial  conditions  (see 
Appendix  A),  the  general  solution  (18)  degenerates  to: 

Tm(t)  = C4>.j(T}eYT  ; y > 0 

In  this  case  any  transverse  initial  displacement  T,.(0)  is 
amplified  according  to  the  multiplier  eYT.  The  real  time  computed 
value  of  9.44  implies  an  amplification  factor  of  nearly  three  orders 
of  magnitude  only  3/4  seconds  after  excitation.  In  practice,  however, 
this  build-up  is  unlikely  owing  to  the  attenuating  effect  of  damping 
and  the  improbable  state  of  tuning  and  initial  conditions  necessary 
for  maximum  y. 

The  probability  of  experiencing  single  round  parametric  resonance 

can  be  significantly  increased  if  the  duration  of  the  ballistic  cycle 

becomes  briefer  than  that  assumed.  Actually,  it  is  the  ratio  of  the 

period  of  the  ballistic  pjlse  - whatever  its  shape  - to  the  funda- 

1 5 

mental  axial  period  which  is  of  importance.  Figure  6 shows  the 

15Harris,  C.  M. , and  Crede,  C.  E.,  Shock  and  Vibration  Handbook, 

Vol.  I - Basic  Theory  and  Measurements,  McGraw-Hill , 1961 , 8-24 . 
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enormous  Influence  of  tnls  ratio  on  the  value  of  y.  For  example,  a 
haversine  pulse  cf  ratio  0,8  will  solicit  an  axial  response  Cd-|(m®*II 
which  Is  about  twenty-five  times  greater  than  the  previous  case 
(where  the  ratio  was  approximately  3.5).  Since  y directly  depends  on 
q-| rox) , an  amp! Iflcatlon  of  three  orders  of  magnitude  is  realized  In 
less  than  seventeen  cycles  of  axial  vibration  (about  3C  milliseconds). 

In  view  of  this  potential  for  large  y it  is  therefore  important  that 
the  design  of  a weapon  be  such  that  axial  vibration  magnitudes  be  kept 
small  - principally  by  creating  Intentional  mismatch  between  the 
fundamental  axial  period  and  the  period  of  the  ballistic  cycle. 

Referring  again  to  figure  6,  a haversine  ballistic  period  should  be  at 
least  twice  the  fundamental  axial  period.  For  ballistic  pulse  shapes 
which  deviate  considerably  from  a haversine,  a response  spectrum 
similar  to  that  of  figure  6 can  be  easily  derived  via  computer. 

Before  moving  on  to  a specific  example  in  multiple  round  resonance, 
some  consideration  should  be  given  to  the  manner  in  which  the  excitation 
differs  from  that  considered  in  the  single  round  situation.  As 
previously  stated  it  is  not  the  free  axial  vibrations  but  rather  the 
quasi-static  responses  of  the  axial  modes  which  serve  as  excitation 
of  transverse  vibrations.  For  example  if  P*(t)  is  a single  haversine 
pulse  of  duration  tq: 

i.e.  P*(t)  = P*/2(l  - cos  2ttt/t  ) 
o o 

then  the  solution  to  (17a)  Is: 
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If  rQ  » 2/ j , the  natural  period,  then 

Mj  = -2P*(  ) / ( jn ) 2 (20) 

that  Is,  the  response  is  quasi-static.  In  all  weapons  of  short  tube 
length  - typically  automatic  small  arms  - the  ballistic  period  is  much 
longer  than  the  fundamental  natural  period  of  axial  vibration  and  the 
assumption,  t0  » 2/j  is  justified.  Thus  the  oscillatory  terms  in 
the  solution  (19)  can  be  neglected  with  little  sacrifice  to  accuracy. 

On  the  other  hand  while  tq  may  be  long  compared  to  the  fundamental 
axial  period,  it  is  very  short  when  compared  to  the  periods  of  the 
first  few  transverse  modes.  It  is  tempting,  therefore,  to  consider 
replacing  the  q.(r)  - as  they  appear  in  (17b)  - by  impulsive  type  terms. 

vJ 

It  will  be  shown  that  such  a replacement  sacrifices  little  ir,  the  way 
of  accuracy  and  leads  to  a very  convenient  consolidation  of  results. 
Multiple  Round  Parametric  Resonance  - Impulse  Excitation 

In  1968,  D.  P.  Krajcinovic  and  G.  Herrmann  published  work55  dealing 
with  the  parametric  resonance  of  bars  subjected  to  repeated  impulsive 
compression.  The  equation  studied  by  the  authors  can  be  obtained  by 
substituting  the  quasi-static  solutions  of  (17a)  into  (17b): 

fM  * *V[1  • P*(t)]TH  * 0 <2>  > 

BKrajcinovic,  P.  P.  and  Herrmann.  G. , Stability  of  Straight  Bars 
Subjected  to  Repeated  Impulsive  Compression,  mIAA  Journal , Oct.  68, 
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where 


and 


'M  = 


fy2  = a%" 

The  load  function  considered  by  Krajcinovlc  and  Herrmann: 

k=oo 

P*(t)  = P,*  + p*  l 6(ex  - k0x>  (22) 

k=-°° 

where  t Is  the  period  of  P*(t).  Physically,  (22)  represents  an 
Infinite  series  of  Impulsively  applied  compressive  loads  superposed 
upon  a steady  load  of  magnitude  P-|*.  6(z)  represents  the  Dirac 

function  and  6 Is  a parameter  having  the  dimensions  of  frequency.  In 
practice  the  authors  did  not  have  to  deal  with  the  full  load  expression 
(22)  but  only  one  cycle  thereof,  for  example  - a single  Impulse  applied 
at  x = t.  This  Is  possible  because  Floquet  theory  (see  Appendix  A) 
enables  the  solution  over  one  period  of  the  excitation  to  be  extended 
Indefinitely  In  time.  Furthermore,  questions  of  stability  can  be 
answered  by  considering  whether  the  motion  grows  or  decays  as  a result 
of  a single  load  application.  Thus  the  load  function  actually  used 
for  analysis  was  equivalent  to: 

P*(t)  = Po*6[0(t  - t)]  (23) 

that  Is,  a single  Impuls ively  applied  load  at  t * r . 

Defining  v^P0*  = u^,  equation  (21)  becomes: 

T + n2[ 1 - u6(6(t  - f ))]T  = 0 (24) 

(The  subscript  M is  implied  throughout.) 
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Stability  Analysis 

Roquet's  theorem'^  guarantees  the  existence  of  particular 
solutions  7(t),  to  equation  (24)  such  that 

T(t  + t)  = pT(t) 


(25) 


where  t is  the  period  of  the  applied  load.  We  are  interested  in 
determining  the  conditions  under  which  p may  be  real  and  have  an 
absolute  value  greater  than  unity  Indicating  a growth  of  the  response 
amplitude  after  one  period  of  the  excitation.  Bolotin^  has  shown 
that  if  two  linearly  independent  solutions  T-j(t)  and  T2(t)  are 
chosen  which  satisfy  the  so-called  'unitary  conditions': 


T-|  (0)  = 1 1,(0)  = 0 

T2(0)  = 0 t2(o)  = 1 


(26) 


then  the  equation  for  p becomes  simply: 

p2  - 2Ap+  1 = 0 

where  .. 

A =1[T,(?+)  + T2(^)] 

The  condition  that  a value  of  p exists  with  |pj  >1,  therefore  amounts 
to: 

| A | > 1 (27) 

It  is  also  apparent  that  the  two  roots  of  the  quadratic  equation  are 
reciprocal  - i.e.,  p,p2  = 1. 


^6Meirovi tch,  L.,  Methods  of  Analytical  Dynamics,  McGraw-Hill,  1970,264. 

^Bolotin,  V.  V.,  The  Dynamic  Stability  of  Elastic  Systems,  Holden-Day, 
1964,  14. 
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Now  two  linearly  independent  solutions  to  equation  (24)  which 
satisfy  the  unitary  initial  conditions  (26)  are: 

T-j  (t)  = cos  frr 
T2(t)  = ^sin  Jh- 


(28) 


. . .T  < T 


In  order  to  use  the  inequality  (27),  the  values  of  Ti(t+)  and 
T2(t+)  are  required.  While  the  functions  themselves  are  continuous 
at  r = t , their  derivatives  are  not  - owing  to  the  application  of  the 
impulse  at  this  instant.  Thus  the  derivative,  T2(t)  cannot  be 
obtained  by  simply  substituting  Into  the  derivative  expression  for  T2. 
The  step  change  in  the  derivative,  however,  can  be  computed  by  a 
direct  integration  of  equation  (24)  over  a short  time  Interval 
containing  the  point  t = t. 


i.e. 


T+e  x+e 

_ * /.  fi2[y6(e(T  - i)-l)]T(r)dT 

T-e  T-e 


Since  T(t)  is  continuous  at  t,  the  term  ft2T  contributes  nothing 
to  the  integral  as  e is  made  vanishingly  small.  Thus 

T+e  n2uT(T) 


T-e 


0 


(e  -*■  0) 


(29) 


The  right  hand  side  of  (29)  derives  from  the  relation: 

/ <S(az )<t>(z)dz  = —■  / <5(z)<f>(-)dz  = ~ <j>(0) 

-oo  | a | .oo  3 1 3 I 

where  <f>(z ) is  a test  function.  Since  9 is  arbitrary  it  can  always  be 
chosen  as  positive  and  the  absolute  signs  omitted  accordingly. 
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Thus,  from  (29): 


r2(T+)  = T2(t)  + T2(t")  = ^ T2(t)  + cos  fix 


-\  = Rj L 


hence 


|A!  = J 


_ fiu 

2cos  fix  +0  sin  fix 


Thus  the  instability  condition  (27)  becomes: 


A | = sin  fix  + cos  fix  | > 


(30) 


In  the  work  performed  by  Krajcinovic  and  Herrmann,  0 is  defined 
as  being  inversely  proportional  to  x (8  * y ) which  is  the  customary 
definition  when  treating  continuous  excitation  functions.  Adhering 
to  this  definition  in  the  case  of  repetitive,  discontinuous  loads, 
however,  leads  to  a load  function  which  is  physically  improbable. 

For  the  load  function  (23): 
i .e. 

the  impulse  of  this  load  under  the  above  definition  of  0 becomes: 


P*(t)  - P *S[6(x  - t)1 
o 


P * 
0 


Po^ 


P0*  / 6[0(t  - r)]dx  = e = 2tt 


(31) 


that  is,  the  impulse  strength  is  seen  to  depend  on  x , the  period  of 

time  between  impulses.  This  is  not  the  physical  situation  of  interest 

in  armament  - and  many  other  - lications  which  imply  a sequence  of 

equally  strong  impulses  independent  of  their  spacing  in  time. 

Fortunately,  these  cases  can  be  handled  merely  by  redefining  0,  i.e., 

let  6 = t » where  xn  represents  the  duration  of  the  ballistic  load 
To  u 

function  as  shown  in  figure  c. 
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Figure  c.  - Periodic  Ballistic  Pressure  Function 


Under  this  definition  of  6,  the  impulse  (31)  becomes: 

P * 

° s p *T  = T 
0 O T0  - *0 


I0  is  defined  as  an  impulse  conveyed  by  a load  of  average  value  PQ* 
and  duration  tq  . In  practice  P*(T)  in  equation  (21)  is  to  be 
replaced  by  the  load  function  Po*<5[9(t-t)]  which  has  the  same  impulse 
as  P*(t).  Thus  Pq*  must  represent  the  nondimensional  time  average 
axial  load  on  the  breech  due  to  ballistic  pressure. 

The  result  of  redefining  9 on  the  instability  condition  (30)  is 
considerable.  Whereas  the  expression  for  A in  (30)  is  purely  a 
function  of  y and  the  ratio  of  the  natural  frequency  to  that  of  the 
excitation,  this  is  not  the  case  when  the  meaning  of  9 is  changed. 
Substituting  9 = 1/tq  into  (30)  the  condition  for  instability  becomes- 


I 


, . . i^llTo  . . . 

1 A J - !~2 — sin  Ur  + cos  Qt|  > 1 


(32) 


Thus  for  the  Mth  equation  (17b)  we  have  p = um  and  U = fiM  as  previously 
defined  so  that  (32)  becomes  finally. 


ig  ' '1  ^ si"  + cos  1 

M^sln  2^+cos  2,31  I 


(33) 


.It 

i-V 

M 1 


It- 

n 

it 

v: 

i: 

S 
! Av 


•fci 

f- 


where  definitions  for  CM  and  w are  inferred.  ue  now  represents  the 
frequency  of  the  impulse  excitation  instead  of  0.  is  seen  to 
depend  on  ld/a  - the  product  of  the  impulse  and  the  slenderness  ratio  - 
and  the  ratio  of  the  unperturbed  natural  frequency  to  that  of  the 
excitation.  It  is  also  apparent  that  in  contrast  to  (30),  AM  is  now 

periodic  in  the  frequency  ratio  implying  equal  stability  criteria  for 

uQm 

all  'zones'  of  instability — = 1/2,  1,  3/2,  ...,  n/2,  ... 

we 

The  Characteristic  Exponent  - Maximum  Value 

The  general  solution  to  Hill's  equation  (17b)  is  given  by  (18). 
This  solution,  having  property  (25),  implies  the  following  definition 
for  y,  the  characteristic,  or  'growth',  exponent: 

LpM 


Y = 


Y real 


For  a given  excitation  frequency  we,  t is  a specified  quantity  and 
maximum  values  of  |p|  determine  maximum  y>  In  calculating  the  partic- 
ular frequency  ratio  for  which  |p|  will  be  maximum,  we  first  let 

ty|/ue  = n/2  + 6 
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n 

*j  j 

ii 


i 

«-■ 

ii.. 


n designates  the  zone  of  Instability  while  6 a local  coordinate  to  be 
varied  In  searching  for  a maximum  value  of  |p|  . Substituting  In  (33) 


/y,  = sin  2tt6  + cos  2tt6](-1  )n 

1^5-  = 0 « (-l)n  [CMI0/a  cos  2-rr6  - sin  2tt6] 

so  that  — tan"^  = 6i,  a value  of  6 for  which  AM  is  extreme. 

2tt  a ” 


From  (26b): 


M 


- 1 


For  extreme  values  of  p: 


3Au  -1/2 

35  = 96“  ^ 4 W - 1 ) ) = 0 


Thus  when  5 = 6-|  , p will  also  be  extreme.  Substituting: 

, , _ . 1 


'ext  1 cos  2tt6i 


± r 


T 


cos22ttS-| 


-1 


hence 


I i i 1 ± 2tt6i  | (34) 

|Plmax  ■ 1 cos  2tt 6 1 1 W 

where  the  sign  is  chosen  to  agree  with  that  of  6-j.  Note  that  |p| 


'max 


is  confirmed  in  (34)  to  be  completely  independent  of  the  zone  number 

Vo 


n.  For  practical  ranges  of 


a 


6-j  is  small  and  can  be  considered 


a linear  function  of  this  parameter, 
i.e.  2ir6i  = Vo/a 
and  (34)  gives:  } p | 


'max 


V( 


a 


or  - tn|p|max  = VQ/a  * 3 degree  of  accuracy. 


(35) 


Figure  7 is  a contour  plot  applicable  to  all  instability  zones. 

The  contours  are  curves  ln|p|  - constant.  The  plot  shows  that  the 

instability  becomes  broader  in  band  and  that  for  larger  values  of  the 

parameter  R = C^I  /:*  the  values  of  Injpj  become  less  sensitive  to 

dR 

variations  in  6,  i.e.,  the  contour  slopes  ^ become  flatter. 

It  still  remains  to  determine  the  conditions  under  which  a 
repetitive  load  function  can  be  justifiably  replaced  in  Hill's 
equation  by  a repetitive  series  of  impul ses.  It  is  expected  that  these 
conditions  will  not  strongly  depend  on  the  particular  shape  of  the 
function  involved  but  rather  its  duration,  t0  and  the  frequency,  ft  , 
of  the  system  on  which  it  acts.  For  the  sake  of  argument,  therefore, 
a sequence  of  rectangular  pulses  ${x}  is  chosen  as  excitation  to  Hill's 


equation: 


f + ft2  (1  - p$(x))f  - 0 


where 


$(t)  = 1,  0 < t£tq 
= 0,  Tq  < T < T 


<5>(x  + t)  = $(x) 

1 8 

The  corresponding  expression  for  A in  this  case  turns  out  to  be  : 


A = (sin  ftxQ  cos  p-|X0 


r P-j  z+ft2 ) 


cos  ftTQ  sin  p-|T0)sin  ftx  + 


/p^+ft2) 

+ (cos  Pt x.  cos  ftxn  + A sin  p,x^  sin  ftx  )cos  ftx 

1 o 0 2np-,  10  0 

where  p-|  = ft /l -y 

*]  0 ~ 

3o<tm,  V.  V.,  The  Dynamic  Stability  of  Elastic  Systems.,  Holden-Day, 
1964,  18. 
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It  Is  easily  verified  that  for  small  values  of  the  quantity  fir0  , 
the  expression  for  A above  reduces  to  that  In  (32).  Expanding  the 
transcendental  terms  containing  the  parameter  fix0  , and  retaining  up 
to  third  order  terms  In  this  parameter,  (36)  becomes: 


A = ftyx0/2  (l+£ryx0/6)  s infix  + cosftx 

Wo  Vo  <*0,  i 

= — — (1  + — 5-)  sin  frr  + cos  nf 


It  Is  therefore  concluded  that  if 

3a 

T°  Wo 

the  approximation  (33)  will  be  good. 


(37) 


EXAMPLE  II  - Ml 39  AUTOMATIC  WEAPON 

The  20mm  Ml 39  at  first  glance  (Fiqure  3)  would  appear  to  be 
vulnerable  to  transverse  vibration  excitation.  It  has  the  appearance 
of  a long  and  slender  tube  - as  compared  with,  say,  a typical  .30 
caliber  gun  such  as  illustrated  in  reference  (9).  Actually  the 
slenderness  ratio  of  the  two  are  nearly  equal  having  values  of 
approximately  103  and  92  respectively.  In  spite  of  this  our  intuition 
is  not  in  error.  As  previously  demonstrated.  It  is  not  the  slenderness 
ratio  alone,  but  its  product  with  the  non-dimensional  ballistic  impulse 
which  determines  the  strength  of  the  exponential  growth  In  parametric 
resonance.  As  It  turns  out,  this  impulse  (I0)  for  the  20mm  round  is 
nearly  three  times  larger  than  that  computed  for  the  .30  cal  weapon. 

A complete  set  of  calculations  follow  - leading  to  the  evaluation 
of  the  parameters  and  I0/a.  The  possible  values  of  the  growth 
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exponent  y are  then  determinable  from  figure  7 . 

A uniform  tube  approximation  is  assumed: 


DQ  = 1.5  in  - tube  outer  diameter 
D.j  = 0.79  in  - tube  inner  diameter 
L = 72  in  - overall  tube  length 
PAV  = 40,000  psi  - time-average  ballistic  pressure 
t0  = .0015  sec  - effective  duration  time  of  ballistic  pressure 
function 

a = 2.02x10®  in/sec  - extension*!  wave  speed 
E = 30x10®  psi  - Young's  Modulus  of  Elasticity 
If  the  ratio  of  the  unperturbed  natural  frequency  to  the 
excitation  frequency  u>e  is  nearly  an  integral  multiple  of  1/2, 
parametric  resonance  is  possible. 


i.e.,  = a6M2  : 


we  2tt/t 


n = 1,2,3... 


or  t = ntT 
afV 


t is  the  time  between  rounds.  In  real  time,  t = xL/a,  corresponding 
to  a firing  rate  of  915  rounds  per  minute  - a reasonable  value  in 
practice.  From  the  previously  assigned  values: 


1 = LvflTT  = 103.3 


- the  slenderness  ratio 


If  we  choose  M = 1,  that  is,  parametric  resonance  of  the  fundamental 
transverse  mode  of  vibration,  then  for  n = 2, 

x = - = 184  and  t = .066  sec 

uB] 


where  we  have  used  the  value  B-j  = 1.875  - corresponding  to  the 
fundamental  mode  of  a cantilevered  beam. 

The  nondlmenslonal  axial  load  (average)*: 

•4 


p * a 

ro 


nPAVD1 2 


4EA 


5.12x10 


The  duration  of  this  load  In  nondlmenslonal  time: 

to  ■ f *0  ■ 

Thus  the  nondlmenslonal  Impulse  Is: 

I = P *t  = 21. 58x1 0“4 
0 0 0 

Recalling  the  definition  of  CM: 

c = y-^!!i_ 

M j (J**)’ 


For  M = 1,  therefore: 

C1  = l 


Ajll 


= .234 


J CJ7r(l  .875 )]2 

The  parameter  n-|T0  must  satisfy  the  smallness  criterion  (38);  i.e., 

filTo  F^t~ 

1 o 

Substituting  values  for  the  parameters: 

.0956  « 57.5 

which  would  appear  to  be  satisfactory. 


♦Note  that  the  real  quantity  of  Interest  Is  the  nondlmenslonal  impulse. 
In.  If  ballistic  curves  are  available  this  quantity  can  be  determined 
directly  by  a simple  integration. 
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Expression  (36),  or  equivalently  figure  7 shows  that  the 
greater f value  for  is  .052.  The  corresponding  growth 

coefficient: 

1nMinax  .052  _ 

'max  ' t .066  ~ *786 

This  value  of  y can  produce  a vibration  amplification  of  over  50 
in  5 seconds  of  repeated  firing.  In  other  words  if  at  any  time  during 
firing  the  deflection  of,  say,  the  muzzle  is  a given  amount  then  5 sec 
later  this  deflection  can  be  over  50  times  greater. 

It  is  not  likely  that  the  exact  state  of  tuning  will  exist  such 
that  the  maximum  amplification  will  be  realized.  Furthermore,  as 
previously  shown,  damping  will  reduce  the  amplification.  On  the  other- 
nand  there  are  other  loads  which  have  not  been  considered  and  which  may 
be  non-mitigating.  Reactions  from  the  moving  projectile  and  the 
travelling  propellant  pressure  as  well  as  the  relative  movements  of 
supports  are  synchronous  with  the  load  considered  in  this  report  and  may 
add  to  its  effect  - possibly  more  than  offsetting  the  reduction  caused 
by  damping.  Again  it  is  mentioned  in  connection  with  figure  7 that 
the  state  of  tuning  leading  to  large  amplifications  need  not  be  as 
precise  when  larger  values  of  the  abcissa  are  encountered. 

The  models  assumed  in  this  report  are  of  course,  somewhat 
idealized.  It  is  always  a question  as  to  when  more  detail  should  be 
’“.eluded.  How  much  is  gained,  for  examole,  by  including  the  geometry 
of  a variable  cross  sectional  area  when  knowledge  of  the  support 


i 
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conditions  Is  always  so  Imprecise?  Further,  It  Is  almost  certain 
that  the  general  character  of  the  results  shown  In  figure  7 will 
prevail  regardless  of  the  detail  Incorporated  In  the  modeling  effort. 
That  Is,  no  matter  how  detailed  the  model,  parametric  resonance  will 
probably  be  confined  to  a narrow  band  of  frequencies  for  the  tensile 
loads  likely  to  be  Induced  on  the  tube  through  firing.  The  frequency 
is  going  to  have  to  be  'just  right'  for  It  to  occur.  But  given  the 
proper  frequency  ratio  our  preceding  analysis  definitely  suggests 
that  ballistic  forces  are  sufficient  to  create  significant  exponential 


growth.  It  Is  doubtful  If  any  higher  degree  of  confidence  in  these 

general  findings  could  be  obtained  through  the  use  of  more  refined 

models.  On  the  other  hand,  the  frequencies  and  even  the  particular 

mode  which  may  resonate  cannot  be  predicted  with  confidence  unless  an 

extremely  detailed  model  Is  employed.  Even  then,  in  cases  where 

parametric  resonance  is  suspect,  field  measurements  should  be  made  to 

determine  precisely  the  frequencies  of  free  vibration  and  excitation. 

As  mentioned  earlier  in  this  report,  there  exists  certain  inconclusive 
9 10 

evidence  * that  dispersion  maxima  observed  in  tests  with  certain 
automatic  weapons  may  be  caused  by  parametric  resonance.  On  the  basis 
of  the  feasibility  demonstrated  herein.  It  Is  advisable  that  such 
experiments  be  repeated  In  a more  tightly  controlled  manner  so  that 
a conclusion  may  finally  be  drawn. 

9AMCP  706-252,  Engineering  Design  Handbook,  Gun  Series,  Gun  Tubes, 
February  1964. 

10Wente,  D.  E.,  SchoenLerger,  R.  L.,  and  Quinn,  B.  E.,  An  Investigation 
of  the  Effort  of  the  National  Frequency  of  Vibration  of  the  Barre 


Upon  the  Dispersions  of  an  Automatic  WeaDon, 


TABLE  I*-  COUPLING  COEFFICIENTS  Ajmn 


f A1n  = -9.83 

■-  ill 

a* 

A121  * 

11.1 

f A112  = 1K1 

A1 22  “ 

-50.1 

[ A113=  7.71 

A123  = 

16.4 

cr> 

11 

A1 24  = 

30.3 

[ A,.  = 1.68 

f lib 

A125  * 

6.93 

1 *116  -8’9 

f' 

A1 26  s 

7.54 

| A,31  ■ 7.7, 

AU1  * 

1.91 

| *132'  16'4 

A142  = 

30.3 

f *133  ■ -,28' 

A1 43  = 

18.4 

[•  *134  ,8'4 

A1 44  = 

-247. 

|i  *136  ' 69'6 

A1 45  = 

18.9 

If  *136  “ ”-3 

4 

A1 46  = 

122 

A|51  = 1,68 

A161 

.819 

A1 52  6,93 

AT62  = 

7.54 

A1 53  = 69,6 

A1 63  = 

11.3 

A1 54  * 19<0 

A1 64  = 

122. 

- A155  = "404* 

A165  = 

19.3 

A1 56  19‘3 

A1 66 

-601. 

*Certain  symmetries  are  evident  throughout  the  Table.  No  attempt  has 
been  made  to  make  use  of  symmetry  properties  in  order  to  shorten  the 
tabulation. 
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I TABLE  1 ■ COUPLING  COEFFICIENTS  A (cont.) 

V 


1 

1 

r 

Agn  = 7.85 

A621  * -24-7 

1' 

i 

A612  = “24*7 

A622  = 101  ’ 

l 

1 . 

L 

f-  - 

A61 3 49‘3 

A623  = -13°- 

L^- 

i 

A614  = "56,3 

A624  = 27°* 

£ 

A615  = 142‘ 

A625  = -878* 

A616  = ~456, 

A626  = 539 • 

£ 

A631  = 49,3 

A641  = _56‘3 

k 

r. 

^632  = '130. 

A642  = 271 • 

f- 

A633  = 296 * 

A643  = "1°78- 

A634  = -1078. 

A644  = 556‘ 

f*  - 

A635  = 552* 

A645  = “432* 

**- 

r 

f. 

u 

A636  = “459‘ 

A646  = 71 6 • 

k 

s 

A651  = 142‘ 

A661  = ~456' 

V 

f* 

■i 

A652  = "878* 

A662  539- 

[? 

A653  552‘ 

A663  = ~459‘ 

If? 

U 

A654  = ”432, 

A664  = 71 6 • 

1 

A655  75°* 

A665  = "632, 

1* 

A656  = "632> 

= g36. 

■ 
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TABLE  II  - EIGENVALUE  ANALYSIS  (NASTRAN)  - MCAAAC , 60MM 
(Trussed  Configuration) 

Mode  Type  Frequency  (hz) 


Recoil  mode 

1.81 

Fundamental 

Transverse  mode 

79.5 

2nd 

Transverse  mode 

100.5 

3rd 

Transverse  mode 

125.6 

4th 

Transverse  mode 

218.7 

5th 

Transverse  mode 

296.2 

6th 

Transverse  mode 

412.7 

7th 

Transverse  mode 

525.9 

Fundamental 

Axial  mode 

580.1 

2nd 

Axial  mode 

1189.2 

3rd 

Axial  mode 

1898.4 

4th 

Axial  mode 

2496.2 
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TABLE  III  - PARAMETRIC  VALUES,  MCAAAC  TUBE 


The  parametric  values  listed  below  are  based  on  the  material  constants 
E,pfor  steel,  a peak  ballistic  pressure  of  75,500  psi,  and  the  NASTRAN- 
predicted  frequencies  for  the  fundamental  axial  vibration  and  the  fifth 
transverse  vibration  modes. 


x/t  = 

£M. 

EA 


6.44xl0'5 

1160 

^(1-  COS  211) 
2 To 

7.03xl0"4 

6.96 
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, .stability  Zones.  Excitation  in  Hill's  Equation  is  a 
Rectangular  Ripple  Function,  (cf  ref  (2),  p.  18,  fig.  4) 
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Figure  5 - Gun  Tube  Model  - An  Axially  Free  Cantilevered  Beam 


Figure  6 - Residual  Response  Spectra,  Haversine  Excitation;  from  Reference  (15) 
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Substituting  Into  A2: 


Af-j  (t+T)  + Bf2(t+T)  * pfAf^t)  + Bf2(t)] 

Using  A3: 

2 

.1  + Ba2jfjCtn  - p[Af,(t)  + Bf2(t)] 

J 1 

Since  f-j  and  f2  are  independent: 


Aa^  + Bu2,  = Ap 
Aa12  + Ba22  = Bp 


A5  only  has  nontrivial  solutions  if  the  determinant  of  the 
coefficients  of  A and  B vanish: 


(A5) 


i .e. 


an  ' p 

a12 


21 

l22 


0 


(A6) 


Substituting  into  either  of  the  equations  A5,  each  root  p. 

J 

gives  a ratio  B/A  such  that  the  solution  will  satisfy  the  condition 
A4. 


d> 

A 0 a21 


(A7) 


In  particular,  if  f^(t)  and  f2(t)  are  chosen  to  satisfy  the  so 
called  'unitary  initial  conditions': 

f](0)  = 1 f2(0)  = 0 

fy  (0)  = 0 f2'  (0)  = 1 


then,  from  A2: 


A = Y(0) 

B = Y’ (0) 
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and  hence  the  ratio  B/A  is  the  ratio  of  initial  conditions  such  that 


A4  Is  satisfied.  Further,  from  A3,  definitions  are  given  to  the  a^: 

a11  ■ fl<T> 

hence 

a21  = f2(T) 

a12  = ^1*  ^ 
a22  = f2*  (T) 


and  [Y1 (0)/Y(C) J , 


Pj  - fl<T> 
f2(T) 


(A8) 


The  quadratic  equation  for  p (A6)  becomes: 


. *★ 
where 


p2  - 2Ap  + B = 0 
A = \ [f]  (T)  + f2'(T)] 


B = f-|  ( T ) f 2*( T ) - fgdJfyd) 


(A9) 


Multiplying  the  first  of  these  equations  by  f2 
and  subtracting: 

(19) 


f f 11 
t1t2 


f 2f 1 M = 0* 


and  the  second  by 


Integrating: 


B = f^tjf^t)  - f.UJf^t)  = const 
The  constant  has  unit  value  in  view  of  the  unitary  initial  conditions 
employed. 


19Boyce,  W.  E.  and  DiPrima,  R.  C.,  Elementary  Differential  Equations 
and  Boundary  Value  Problems,  John  Wiley,  1966,  89.  : 

*Note  that  this  equation  can  be  written  W'  = 0,  where  W19  is  the 
Wronskian  of  f-j  (t)  and  f2(t).  "'2  u + 

**If  any  of  the  fj  have  jump  discontinuities  at  t = T,  then  A = A(t  ) 
is  implied. 


'j2 


ut  \ V**‘t  «*■  > ««.»• 


Thus  equation  A6  becomes: 

p2  - 2Ap  + 1 = 0 


(A10) 


whereupon  _____ 

p = A ± /a  2 ' . l 

Note  that  the  two  roots  of  this  equation  are  reciprocal,  i.e.  p^  = 1 

Thus  corresponding  to  real  roots,  A4  represents  a solution  which 

grows  or  diminishes  (according  to  the  choice  of  p)  following  each 

period  of  the  excitation,  T.  The  ratio  of  initial  conditions  which 

will  lead  to  either  of  these  solutions  alone  is  given  by  A8. 

EXAMPLE  - Delta  Function  Excitation 

Let  00 

$(t)  = l S(t-kT) 
k=l 

From  A10,  the  roots  p.  depend  entirely  on  the  quantity 

J 


A ■ \ [f,(T)  + f2'(T)] 

where  f^(t)  and  f2 (t)  are  linearly  independent  solutions  of  Hill's 
Equation  (A1 ) satisfying  unitary  initial  conditions.  These  can  be 
chosen  as: 

f^(t)  = cos  Lt 

f2(t)  = ?2sinnt  (A11) 

Thus  f-j  (T)  = cos 

The  expression  obtained  by  a formal  differentiation  of  f2(t) 
above,  is  valid  up  to  but  not  including  time  t = T since  the 
application  of  the  first  'impulse'  will  cause  f2‘  (T)  to  be  discontinuous. 

This  jump  discontinuity  can  be  evaluated  by  an  integration  of  Hill's 

✓ 

Equation  (Al): 
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Thus 


r 

£ 

r 

; * . 
t-i 

l , 


and 


T+At 

U (T)  - f ' (T”)  = 11m  /Q2[e6(t-T)-1 ]f9(t)dt  = eQ2f9(T) 
At-MD  T-At  £ 


f2‘  (T)  = efisinnT  + cosftT 


A = — sinfiT  + cosQT 
2 


whereupon  the  roots  p^,  p2  may  be  evaluated. 

For  any  given  finite  value  of  e,  intervals  of  the  parameter  JlT 
exist  such  that  these  roots  will  be  real  and  since  their  product  must 
be  unity,  the  following  possibilities  arise: 

(0  p1  = p2  = ± 1 
(ii)  p1  > 1,  P2  < 1 

Case  (i)  represent  borderline  situations,  that  is,  the  solutions  A4 
neither  grow  or  decay  with  each  period  T of  the  excitation.  Case  (ii), 
on  the  other  hand,  produces  one  solution  of  th;»  form  A4  which  grows  In 
amplitude  by  a factor  p1  > 1 and  another  which  decays  by  the  factor 
P2  < 1,  each  period  of  the  excitation.  The  initial  conditions  which 
produce  either  of  these  solutions  must  only  satisfy  A8, 

= = [Pi  - «s!JT]!!/s1n!)T 

Figure  Al(a,b)  shows  the  solution  when  this  ratio  is  enforced.  In 
general,  however,  the  solution  for  arbitrary  initial  conditions  is  a 
linear  combination  of  these  special  solutions  of  pure  growth  and  decay. 
Figure  Al{c)  shows  one  such  solution  exhibiting  early  decay  even- 
tually to  be  overpowered  by  the  growing  solution. 
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; j 
1 i 


The  solutions  of  pure  growth  or  pure  decay  can  be  represented  as: 


where 

so  that  indeed. 


i In  p.. 
fj(t)  - ^j(t)eT  J 


*j(t+T)  = *j(t) 


(f  + Din  Pi 
f j ( t+T ) - H»j(t)e  T J 


pjfj 


as  required.  In  general  p is  complex 
hence 

In  p.  s ln|p.|  + i arg(p.) 

J J J 

Thus  $■  ln|p.  I 

fj(t)  = 4>j  (t)e'  j 

where  it  arg  p^ 

4>j(t)  = t|»j(t)e'T  J 

When  the  p^.  are  real,  4>j(t)  = ipj (t ) and  are  periodic 

the  general  solution  to  Hill's  equation  may  be  written 

7 1 n | p-j  | ~ ln|p2| 

Y(t)  = C1<}>1  (t)eT  + C2q>2(t)eT 

or,  enforcing  the  relation  P-|P2  = 1: 


Y(t)  = C-j  <f>-j  (t)e 


f lnIPil  + 


-t 


C2<t>2(t)e 


Tln|pl 


^Bolotin,  V.  V. 
1964,  14. 


The  Dynamic  Stability  of  Elastic  System; 
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(A12) 

(t) 

(A13) 

in  T.  Thus 


, Holden-Day, 


FIGURE  A1 


- SOLUTIONS  TO  HILL'S  EQUATION  - IMPULSE  EXCITATION  OF 
STRENGTH  0.1  APPLIED  TO  A SYSTEM  OF  UNIT  CIRCULAR 
FREQUENCY.  INITIAL  CONDITIONS  INDICATED  IN  THE 
FOLLOWING  CAPTIONS  (Al(a),  Al(b),  Al(c)  ) ARE  PREDI- 
CATED ON  AN  ASSUMED  INITIAL  ENERGY  PER  UNIT  MASS  OF  1/2. 


y(0)  = -.6892250 


00*0  ovo 
lN3W33tndSI0 


;>9 


v: 


Fi  qi-re  A1  (c)  - Gerv 
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